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Abstract. We discuss the problem of finding an upper bound for the number 
of equilibrium points of a potential of several fixed point charges in R n . This 
question goes back to J. C. Maxwell I1UI and M. Morse 1121 . Using fewnomial 
theory we show that for a given number of charges there exists an upper bound 
independent of the dimension, and show it to be at most 12 for three charges. 
We conjecture an exact upper bound for a given configuration of nonnegative 
charges in terms of its Voronoi diagram, and prove it asymptotically. 
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1. Introduction 

Consider a configuration of I = ji + v fixed point charges in W , n > 3 consisting 
of [i positive charges with the values Ci> • • • j Cad an d v negative charges with the 
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values ■ ■ ■ >0- They create an electrostatic field whose potential equals 



where r, is the distance between the i-th charge and the point x = (xi , . . . ,x n ) G R™ 
which we assume different from the locations of the charges. Below we consider the 
problem of finding effective upper bounds on the number of critical points 
of V(x), i.e. the number of points of equilibrium of the electrostatic force. In 
what follows we mostly assume that considered configurations of charges have only 
nondegenerate critical points. This guarantees that the number of critical points is 
finite. Such configurations of charges and potentials will be called nondegenerate. 
Surprisingly little is known about this whole topic and the references are very scarce. 

In the case of R 3 one of the few known results obtained by direct application of 
Morse theory to V{x) is as follows, see ^2]) Theorem 32.1 and Theorem 6. 

Theorem 1.1 (Morse-Kiang). Assume that the total charge 2i=i Cj * n fZHP * s 
negative (resp. positive). Let mi be the number of the critical points of index 1 of 
V , and ni2 be the number of the critical points of index 2 ofV. Then mi > /i (resp. 
rri2 > fJ> — 1) and mi > v — 1 (resp. mi >v). Additionally, mi — m2 = v — ji— 1. 

Note that the potential V(x) has no (local) maxima or minima due to its har- 
monicity. 

Remark. The remaining (more difficult) case X)f=i + X^- =AI +i Cj = is treated 
in 0. 

Remark. The above theorem has a generalization to any R™, n > 3 with mi being 
the number of the critical points of index 1 and rr»2 being the number of the critical 
points of index n — 1. 

Definition 1.2. Configurations of charges with all nondegenerate critical points 
and mi + m2 = /i + v + 1 are called minimal, see |12| . p. 292. 

Remark. Minimal configurations occur if one, for example, places all charges of the 
same sign on a straight line. On the other hand, it is easy to construct generic 
nonminimal configurations of charges, see |12) . 

Remark. The major difficulty of this problem is that the lower bound on the 
number of critical points of V n given by Morse theory is known to be not exact. 
Therefore, since we are interested in an effective upper bound, the Morse theory 
arguments do not provide an answer. 

The question about the maximum (if it exists) of the number of points of equi- 
librium of a nondegenerate configuration of charges in R 3 was posed in ^2]i P- 293. 
In fact, J. C. Maxwell in ^U], section 113 made an explicit claim answering exactly 
this question. 

Conjecture 1.3 see also 2]below). The total number of points of equilibrium 
(all assumed nondegenerate) of any configuration with I charges in R 3 never exceeds 



Remark. In particular, there are at most 4 points of equilibrium for any configura- 
tion of 3 point charges according to Maxwell, see Figure 



(1.1) 




(I -If. 
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Before formulating our results and conjectures let us first generalize the set-up. 
In the notation of Theorem 11.11 consider the family of potentials depending on a 
parameter a > and given by 

M ^Kf + - + l) + te + "- + f)' 

where pi — rf, i — 1, (The choice of p^'s instead of r^s is motivated by 

convenience of algebraic manipulations.) 

Notation 1.4. Denote by Ni(n,a) the maximal number of the critical points of the 
potential M.fy) where the maximum is taken over all nondegenerate configurations 
with I variable point charges, i. e. over all possible values and locations of I point 
charges forming a nondegenerate configuration. 

Our first result is the following uniform (i.e. independent on n and a) upper 
bound. 

Theorem 1.5. a) For any a > and any positive integer n one has 
(1.3) Ni(n,a) <4 l \3l) 21 . 

b ) For I — 3 one has a significantly improved upper bound 

N 3 (n,a) < 12. 

Remark. Note that the right-hand side of the formula (|1.3|l gives even for I = 3 the 
horrible upper bound 139, 314, 069, 504. On the other hand, computer experiments 
suggest that Maxwell was right and that for any three charges there are at most 4 
(and not 12) critical points of the potential (|1.2fl . see Figure^ 




Figure 1. Configurations with two and with four critical points. 

Remark. Figure ^ shows the level curves of the restrictions of the potential of three 
positive charges to the plane they span in two essentially different cases (conjec- 
turally, the only ones). The graph on the left has 3 saddles and 1 local minimum 
and the graph on the right has just 2 saddle points. 
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1.1. Voronoi diagrams and the main conjecture. Theorem 11.61 below deter- 
mines the number of critical points of the function V a for large a in terms of the 
combinatorial properties of the configuration of the charges. To describe it we need 
to introduce several notions. 

Notation. By a (classical) Voronoi diagram 1 of a configuration of pairwise 
distinct points (called sites) in the Euclidean space R n we understand the partition 
of K" into convex cells according to the distance to the nearest site, see e.g. [2] and 

El 

A Voronoi cell S of the Voronoi diagram consists of all points having exactly 
the same set of nearest sites. The set of all nearest sites of a given Voronoi cell S 
is denoted by J\fS(S). One can see that each Voronoi cell is a interior of a convex 
polyhedron, probably of positive codimension. This is a slight generalization of 
traditional terminology, which considers the Voronoi cells of the highest dimension 
only. 

A Voronoi cell of the Voronoi diagram of a configuration of sites is called effec- 
tive if it intersects the convex hull of J\fS(S). 

If we have an additional afhne subspace L C K™ we call a Voronoi cell S of the 
Voronoi diagram of a configuration of charges in E™ effective with respect to L 
if S intersects the convex hull of the orthogonal projection of AfS(S) onto L. 

A configuration of points is called generic if any Voronoi cell S of its Voronoi 
diagram of any codimension k has exactly k + 1 nearest cites and does not intersect 
the boundary of the convex hull of AfS(S). 

A subspace L intersects a Voronoi diagram generically if it intersects all its 
Voronoi cells transversally, any Voronoi cell S of codimension k intersecting L has 
exactly k + 1 nearest sites, and S does not intersect the boundary of the convex 
hull of the orthogonal projection of JVS(S) onto L. 

The combinatorial complexity (resp. effective combinatorial complex- 
ity) of a given configuration of points is the total number of cells (resp. effective 
cells) of all dimensions in its Voronoi diagram. 

Example. Voronoi diagram of three non-collinear points A,B,C on the plane con- 
sists of seven Voronoi cells: 

(1) three two-dimensional cells Sa,Sb, Sc with MS(S) consisting of one point, 

(2) three one-dimensional cells Sab, Sac, Sbc with MS(S) consisting of two 
points. For example, Sab is a part of a perpendicular bisector of the 
segment [A, B\. 

(3) one zero-dimensional cell Sabc with AfS(S) consisting of all three points. 
This is a point equidistant from all three points. 



The first known application of Voronoi diagrams can be traced back to Aristotle's De Caelo 
where Aristotle asked how a dog faced with the choice of two equally tempting meals could ratio- 
nally choose between the two. These ideas were later developed by the known French philosopher 
and physicist Jean Buridan (1300-1356) who sowed the seeds of religious scepticism in Europe. 
Buridan allowed that the will could delay the choice in order to more fully assess the possible 
outcomes of the choice. Later writers satirized this view in terms of an ass who, confronted by 
two equally desirable and accessible bales of hay, must necessarily starve while pondering a deci- 
sion. Apparently the Roman Catholic Church found unrecoverable errors in Buridan's arguments 
since about hundred and twenty years after his death a posthumous campaign by Okhamists suc- 
ceeded in having Buridan's writings placed on the Index Librorum Prohibitorum (List of Forbidden 
Books) from 1474-1481. 
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There are two types of generic configurations, see Figure First type is of an 
acute triangle AABC and then all Voronoi cells are effective. Second type is of an 
obtuse triangle AABC, and then (for the obtuse angle A) the Voronoi cells Sbc 
and Sabc are not effective. 

The case of the equilateral triangle AABC is non-generic: the cell Sabc, though 
effective, lies on the boundary of the triangle. 

The following result motivates our main conjecture If .71 below. 

Theorem 1.6. a) For any generic configuration of point charges of the same 
sign there exists ao > such that for any a > ao the critical points of 
the potential V a (x) are in one-to-one correspondence with effective cells of 
positive codimension in the Voronoi diagram of the considered configuration. 
The Morse index of each critical point coincides with the dimension of the 
corresponding Voronoi cell. 
b) Suppose that an affine subspace L intersects generically the Voronoi diagram 
of a given configuration of point charges of the same sign. 

Then there exists ao > ( depending on the configuration and L ) such 
that for any a > ao the critical points of the restriction of the potential 
V a (x) to L are in one-to-one correspondence with effective w.r.t. L cells of 
positive codimension in the Voronoi diagram of the considered configuration. 
The Morse index of each critical point coincides with the dimension of the 
intersection of the corresponding Voronoi cell with L. 

Finally, our computer experiments in one- and two-dimensional cases led us to 
the following optimistic 

Conjecture 1.7. a) For any generic configuration of point charges of the 
same sign and any a > h one has 

(1-4) < < f , 

where a? a is the number of the critical points of index j of the potential 
V a (x) and \ 3 is the number of all effective Voronoi cells of dimension j in 
the Voronoi diagram of the considered configuration. 
b) For any affine subspace L generically intersecting the Voronoi diagram of 
a given configuration of point charges of the same sign one has 

(1-5) <£ < Hi, 

where a J a L is the number of the critical points of index j of the poten- 
tial V a (x) restricted to L and jj^ is the number of all Voronoi cells with 
dim(S fl L) = j effective w.r.t L in the Voronoi diagram of the considered 
configuration. 

We will refer to the inequality (|f .4(1 resp. fl .5(1 as Maxwell resp. relative 
Maxwell inequality. 

Remark. Theorem If .61 and Conjecture If .71 were inspired by two observations. On 
one hand, one can compute the limit of a properly normalized potential V a (x) when 
a — > oo. Namely, one can easily show that 

lim V a " (x) = Voo(x) = min p l (x). 
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This limiting function is only piecewise smooth. However, one can still define critical 
points of Voo(x) and their Morse indices. Moreover, it turns out that for generic 
configurations every critical point of Voo(x) lies on a separate effective cell of the 
Voronoi diagram whose dimension equals the Morse index of that critical point, 
see 12.4.31 Theorem If .61 above claims that for sufficiently large a the situation is 
the same, except that the critical point does not lie exactly on the corresponding 
Voronoi cell (in fact, it lies on 0(a~ 1 ) distance from this Voronoi cell, see Lemma 
12.291 and I2.34[) . On the other hand, computer experiments show that the largest 
number of critical points (if one fixes the positions and values of charges) occurs 
when a — » 00. 

Even the special case of the conjecture If .71 when L is one-dimensional is of in- 
terest and still open. Its slightly stronger version supported by extensive numerical 
evidence can be reformulated as follows. 

Conjecture 1.8. Consider an l-tuple of points (x\, y\), . . . , (xi,yi) in R 2 . Then 
for any values of charges (£1, . . . , Q) the function V*(x) in (one real) variable x 
given by 



Remark. In the simplest possible case a = 1 conjecture 1 1. 81 is equivalent to showing 
that real polynomials of degree (4/ — 3) of a certain form have at most (21 — 1) real 
zeros. 

1.1.1. Complexity of Voronoi diagram and Maxwell's conjecture. In the classical 
planar case one can show that the total number of cells of positive codimension of 
the Voronoi diagram of any / sites on the plane is at most 5/ — f 1 and this bound 
is exact. 

Since (/ — l) 2 is larger than the conjectural exact upper bound 5/ — 11 for all 
I > 5 and coincides with 5/ — 11 for I = 3, 4, we conclude that Con i ecture 1 1 . 71 implies 
a stronger form of Maxwell's conjecture for any I positive charges on the plane and 
any a>\. 

For n > 2 the worst-case complexity T(l,n) of the classical Voronoi diagram 
of an Z-tuple of points in R n is 8(^ n / 2+1 ]), see 0. Namely, there exist positive 
constants A < B such that Al^ n / 2+ ^ < T(l,n) < Bl^ n / 2+1 l Moreover, the Upper 
Bound Conjecture of the convex polytopes theory proved by McMullcn implies that 
the number of Voronoi cells of dimension k of a Voronoi diagram of / charges in R™ 
does not exceed the number of (n — fc)-dimensional faces in the (n + l)-dimensional 
cyclic polytope with I vertices, see This bound is exact, i.e. is achieved for 

some configurations, see |15j . 

In R 3 this means that the number of O-dimensional Voronoi cells of the Voronoi 
diagram of I points is at most ; ('~ 3 ) , the number of 1-dimensional Voronoi cells is 

at most 1(1 — 3), and the number of 2-dimensional Voronoi cells is at most . 

We were unable to find a similar result about the number of effective cells of 
Voronoi diagram. However, already for a regular tetrahedron the number of effective 
cells is 11, which is greater than the Maxwell's bound 9. Thus a stronger version 
of Maxwell's conjecture in R 3 fails: the number of critical points of V a could be 
bigger than (I — l) 2 for a sufficiently large. 



(1.6) 




c 



has at most (21 — 1) real critical points, assuming a > h. 
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However Maxwell's original conjecture miraculously agrees with the Maxwell 
inequalities (|1.4(l and we obtain the following conditional statement. 

Theorem 1.9. Conjecture \1.7\ implies the validity of the original Maxwell's conjec- 
ture for any configuration of positive charges in R 3 in the standard 3 -dimensional 
Newton potential, i.e. a = i. 

Existing literature and acknowledgements. Logarithmic potentials in R 3 similar to 
(11.2(1 (i.e. the case of the electrostatic force proportional to the inverse of the 
distance) were studied in a number of papers of J. L. Walsh, see {F^ and references 
therein. In this case it is possible to generalize the classical Gauss-Lucas theorem 
and some results on Jensen's circles for polynomials in one complex variable to real 
vector spaces of higher dimension. 

Critical points of a logarithmic potentials in R 2 — C are zeros of some univariate 
polynomial of degree < I — 1, so the upper bound for the number of critical points 
is Z — 1, see UJ. 

Some interesting examples of electrostatic potentials whose critical points form 
curves were considered in 4 . The question whether degenerate electrostatic poten- 
tial defined by a finite number of charges can have an analytic arc of critical points 
was stated in ^2], P-294. Finally, the results about instability of critical points 
for more general potentials and dynamical systems are obtained in 0. Instability 
in our context follows from subharmonicity of the considered potential and was 
already mentioned in ^U], section 116 under the name Earnshaw's theorem. 

The structure of the paper is as follows. In fJ21we prove the above results and 
present the computer evidence for our main conjecture, [^contains further remarks 
and open problems related to the topic. Finally, in 2] we reproduce the original 
section 113 of ^U] where Maxwell presents the arguments of Morse theory (devel- 
oped at least 50 years later), and names the ranks of the 1st and the 2nd homology 
groups of domains in R 3 in the language of (apparently existing) topology of 1870's 
to formulate his claim. 

The authors are sincerely grateful to A. Eremenko, A. Fryntov, D. Khavinson, 
H. Shapiro, M. Shapiro and A. Vainshtein for valuable discussions and references. 

2. Proofs 

We start this section with a discussion of the nondegeneracy requirement and 
the (co)dimension of the afhne span of a configuration of point charges. 

2.1. Relation between number of charges and dimension. Consider a non- 
degenerate configuration of I = fi + v point charges in M", and let L C R™ be the 
afhne subspace spanned by the points where the charges are located. Evidently, 
dimL <l-l. 

Theorem 2.1. If all critical points of the potential V a are isolated, then either all 
critical points belong to L or n < I — 1. 

Let us first show that it is enough to consider the cases n < I only. 

Lemma 2.2. If a configuration of charges in R n has only isolated critical points 
then either all its critical points belong to L or L is a (real) hyperplane in R". 

Proof. Indeed, assume that there is a critical point outside L and codimL > 1. 
Then the whole orbit of this point under the action of the group of rotations of R™ 
preserving L consists of critical points (since this action preserves the potential) . □ 
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To complete the proof one has to exclude the case n = I. We show that if 
dimL = 1 — 1 then all critical points of the potential are in L. 

Lemma 2.3. If one can find a hyperplane H in L separating positive charges 
from the negative ones, then the potential of the configuration has no critical points 
outside L. 

Proof. Indeed, let x $ L be any point outside L, and let H x be any hyperplane 
containing both x and H and transversal to L. Let n be a vector normal to H x at 
x. The the signs of scalar products of the gradients of the potentials of each charge 
with n are the same, so x cannot be an equilibrium point. □ 

Corollary 2.4. Potential of any configuration of positive charges has no critical 
points outside L. 

Corollary 2.5. Any configuration with I point charges such that dimL = 1 — 1 has 
no critical points of potential outside L. 

Proof. These points should form a non-degenerate simplex, and any subset of ver- 
tices of a simplex can be separated from the rest of the vertices by a hyperplane, 
so the claim follow from the previous Lemma. □ 

Remark. As one can see from the proof, the set of critical points of a configuration 
is a union of spheres with centers in L and of dimension equal to codimL. As n 
grows, the change of the dimension of spheres is the only parameter that changes, 
so the case codimL = 1 is the most general one. 

We conclude that in any case it is enough to consider the case n < I — 1. 

2.2. Proof of Theorem ll.51 a. The proof is an application of the theory of fewno- 
mials developed by A. G. Khovanskii in 0. A serious drawback of this theory is 
that the obtained estimates, though effective, are usually highly excessive. Apply- 
ing the methods, rather than the results of this theory one might get a much better 
estimate which we illustrate while proving Part b) of Theorem II .51 
We start with the following result from jSJ §1.2]. 

Theorem 2.6 (Khovanskii). Consider a system of m quasipolynomial equations 

Pl(u,w(u)) = ... = P rn (u,w(u)) = 0, U = (ui, ...,u m ), 

where each Pi is a real polynomial of degree di in (m+k) variables (ui, u m ,wx, ...,Wk) 
and 

wj = exp(aj, u), a ■ = (a], a" 1 ) e lR m , j = l,...,k. 
Then the number of real isolated solutions of this system does not exceed 

d 1 ---d m (d 1 + ---+d m + l) k 2 k ^/ 2 . 

The estimate of Theorem II .51 a will follow from a presentation of the critical 
points of a configuration of point charges as solutions of an appropriate system of 
quasipolynomial equations, see below. 
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2.2.1. Constructing a quasipolynomial system. Consider a configuration with I point 
charges in R n . Denote by x = (xi, ...,x n ) the coordinates of a critical point and 
denote by (c\, c l n ), i = the coordinates of the i-th charge. We assume 

that the 1-st charge is placed at the origin, i.e. that c\ = • • ■ = c\ = 0. 

The first I equations of our system define the indeterminates p = (p±, . . . , pi) as 
the squares of distances between the variable point x and the charges. They can 
be rewritten as 

(2.1) P 1 (x,p) = ... = P l (x,p) = 0, 
where 

n n 

(2.2) Pi(x,p) =J2 x2 j-Pi' Pi{x,p)=pi-p i + ^c){2x j +c i j ), i = 2,...,l. 

The second group of equations expresses the fact that the point x = (x\, x n ) 
is the critical point of the potential V a (x) =J2 C,iP~ a ■ Namely, 

d 1 d 1 

— VaOr) = ^2(iQ^-Pi a = -2a^2(m(xj - c)) = Pi +j (x,v), j = 1, ...,n 

3 i—l 3 2—1 

where we denote 

(2.3) v i =p~ a ~ 1 , i=l,...,l and v = (v%, . . . ,vi). 
Introducing variables Si = log pi we get the system: 

Pi(x,s,p,v) = ... = P n+ i(x,s,p,v) = 0, 

of (I + n) quasipolynomial equations in (I + n) variables (x, s, p, v) with 

Pi — exp(sj), Vi = exp(-(a + l)s,/2), i=l,...,l. 

This system has the type described in Theorem 12.61 with m = n + /, k = 21, 
degPi = degP i+ i = ... = degP n+k = 2, and degF 2 = •■■ = degF; = 1. By 
Proposition 12 . II one has n < I — 1 which implies the required estimate: 

Ni(n,a) < N t (l - l,a) < ^t^l 21 = 4^(3Z) 2Z . 
For example, for I = 3 one gets N 3 (n, a) < 139, 314, 069, 504. 

2.3. Proof of Theorem ll.51 b. As we mentioned above, one can do much better by 
applying the fewnomials method rather than results, and here we demonstrate this 
in the case of three charges. By Proposition 12 . II we can restrict our consideration 
to the case of R i_1 = M 2 . We also use coordinates (x,y) instead of (xi 7 x 2 ). 

The scheme of this rather long proof is as follows. We make a change of variables, 
passing from (x,y) to new variables (/, g). In the new coordinates the equilibrium 
points coincide with the intersection points of two explicitly written planar curves 
71 and 72 in the positive quadrant R 2 ^ of the real plane. Both curves are separating 
solutions of Pfaffian forms. The Rolle-Khovanskii theorem applied twice produces 
two real polynomials R and Q such that the required upper bound can be given 
in terms of the number of their common zeros in . The latter is bounded from 
above by the Bernstein-Kushnircnko bound minus the number of common roots of 
R and Q lying outside R+. 
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2.3.1. Changing variables and getting system of equations. To emphasize that the 
methods given below can be generalized we make the change of variables in the 
situation of I charges in 1" (we assume, as before, that n < I — 1). This, as a 
byproduct, produces another proof of Theorem II . 51 with a somewhat better upper 
bound. 

As above, we assume that the charges Q are located at (c\, ■ ■■c\ l ), i = 1, I, and 
consider the potential 

l n 

V a (xi, ...,x n ) = '^2dp~ a , where pi = / ](xj - c)) 2 , i = l,...,l. 

i=l j=l 

The system of equations defining the critical points of V a (x) is 

dV a (x) . dV a (x) ^ -a-xt i\ 

fl , = 0, j = l,...,n, where w = -2a^Q Pi (»j ~ <$)• 

Introducing hi = p.^"" 1 one can solve each equation of this system and express 
Xj in terms of he 

i i 

(2.4) Xj = — , where a = }^(ihj, (Tj = }^ Qhjc) 

a i=l i=l 

are homogeneous linear functions of hi . The equilibrium points correspond to the 
solutions of the following system of equations obtained from the definition of hi by 
substitution of Ojja instead of Xj: 

j_ t. " 

(2.5) Q+1 = ^, where & = ^(a, - c)af , i = 1, ...,l. 

G 3=1 

This system has following remarkable properties: 

Proposition 2.7. a) Any solution of a = = is also a zero of all £i 's; 
b ) each £j is a strictly positive real quadratic polynomial independent of hi . 

Proof. Indeed, & - & - c})(2 CTj - a(cj + c))). 

The second statement is evident except the independence on hi, which is proved 
by direct computation. □ 

Remark. Note that the above system can be represented as a system of quasipoly- 
nomials as in Theorem 12.61 Namely, the equations in i|2.5|) are polynomials in 
hi, hl^ a+1 ^ . Introducing Si — log fi one can apply Theorem l2.6l After several small 
tricks - dehomogenization of the system, introduction of a new variable z = £i and 
noting that the expression for £i — £j becomes then linear - we obtain an upper 
bound 2 ■ A 1 (21 + 3) 21 on the number of equilibrium points of a system of I charges. 
For I > 3 this bound is somewhat better than the bound fOl 

Now, let us use the previous construction for I — 3 and n — 2. Without loss of 
generality we can assume that the three charges with the values £i, £2, 1 are located 
at (0,0), (1,0) and (a, b) respectively. 

Expressions l|2.4|l are homogeneous in hj, so we introduce the nonhomogeneous 
variables / and g as follows 

, / \ a+l , / \ a+1 

(2.6) f=^=( P A and g^-f* 



hi \p2 J hi \p 3 
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Then equations (|2.4|l become: 

(2 7) ~- a 9 + (zf .._ b 9 



C1+C2/ + 5' C1+C2/ + 5' 

The system (|2.5|1 reduces to the following two equations describing two curves 
71 and 72 in the positive quadrant of the (/, g)-plane: 

(2.8) 7i = {/ 1/(Q+1 W = l}, 7 2 = {ff- 1/(tt+1) 6 = r 1/(Q+1) 6}. 
Here 

6 = («.9 + C 2 /) 2 + foV, 

& = (( a -i)<?-Ci) 2 + &v, 

6 = ((a-l)C 2 / + ad) 2 + 6 2 (C 2 / + Ci) 2 . 
The following facts about £j follow from the Proposition 12. 71 

Proposition 2.8. (1) £2 depends only on g, and £3 depends only on f; 

(2) £2,^3 are strictly positive quadratic polynomials; 

(3) £1 is a positive definite homogeneous quadratic form; 

(4) £1 = £2 = £3 = Ziawe £iajo complex solutions. 

The goal of all subsequent computations is to give an upper bound on the number 
N of the points of intersection of 71 and 72 in M%_. We are able to obtain the 
following estimate proved below. 

Proposition 2.9. The number of intersection points 0/71 and 72 lying in is 
at most 12. 

Note, that any intersection point of 71 and 72 lying in the positive quadrant 
= {/ > 0,g > 0} corresponds, via (|2.7|) . to a unique critical point of V a (x,y), 
so the estimate of II .51 b immediately follows. 

2.3.2. Rolle-Khovanskii theorem. Before we move further let us recall the IR 2 -version 
of a generalization of Rolle's theorem due to Khovanskii. 

Suppose that we are given a smooth differential 1-form u defined in a domain 
D C M 2 . Let 7 C D be a (not necessarily connected) one-dimensional integral 
submanifold of u. 

Definition 2.10. We say that 7 is a separating solution of w if 

a) 7 is the boundary of some (not necessarily connected) domain U; 

b) the coorientations of 7 defined by lu and by U coincide (i.e. lu is positive 
on the outer normal to the boundary ofU). 

Let 71, 72 be two separating solutions of two 1-forms u>\ and U2 resp. 

Theorem 2.11 (see [SJ). 

tt(7i,72) < |l(7i) + KTijTa), 

where $(71,72) is the number of intersection points of 71 and ^2, tt(7i) *s the number 
of non- compact components 0/71 and t»(7i, 72) is the number of the points of contact 
of 71 and L02, i.e. the number of points 0/71 sucft i/iai £^2(71) = 0. (One can also 
characterize the latter points as the intersection points of 71 with an algebraic set 
{ux/X^ = 0}.) 
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2.3.3. First application of Rolle-Khovanskii theorem. We apply Theorem 12.111 to 
the curves 71,72 defined in 1)2.8(1 . These curves are integral curves in of the 
one-forms r/\ and 772 respectively, where 

(2.9) m = r *+&°-!!L t 

1 , &\ M , ( 1 & 



(2 - 10) ^ = h^ + gj*H<^6J*- 

These forms are logarithmic differentials of the functions defining the curves: if 
we denote F = (//.g)~ 1/(Q+1) 6^\ thcn 12 = {F = 1} and 772 = dlogF. Similarly, 
th = dlogG, where G = /VC"* 1 )^ 1 . 

In what follows, we assume that 1 is a regular value of F and G. This can be 
always achieved by a small perturbation of parameters and is enough for the proof 
of Theorem ll.5l b by upper-continuity of the number of the non-degenerate critical 
points. 

Lemma 2.12. The curves 71 and 72 are separating leaves of the polynomial forms 
771 and 772. 

Proof. Indeed, 72 is a level curve of the function F = (//s) _1/(a+1) &6f which 
is a smooth function on R+. Thus, 72 coincides with the boundary of the domain 
d{F < 1}. Therefore the value of 772 = d(\ogF) on the outer normal to {F < 1} is 
non-negative, and is everywhere positive since 1 is not a critical value of F. This 
means that the coorientations of 72 as the boundary of {F < 1} and as defined by 
the polynomial form 772 coincide. 

Similar arguments hold for 71, and we conclude that 71 and 72 are separating 
leaves of the forms 771 and 772. □ 

This enables application of Theorem 12 . 1 II to the pair (71,72) and we obtain the 
following estimate: 

Proposition 2.13. 

N < Ni + N 2 , 

where N is the number of points in the intersection 71 072 nl^_, iVi is the number 
of the noncompact components of 72 in R+ and N 2 is the number of points of 
intersection 0/72 with the set T = {7/1 A 772 = 0} in 

Lemma 2.14. Ni = 2. 

Proof. Asymptotically the equation / _1 '' a+1 ^3 = 1 ^ ( - Q+1 ^2 has four solutions: 
g ~ const •/ as / — > or 00, g ~ const ■f^ 1 ~ 2a as / — > 00, and / ~ const -(? _1_2q 
as g — y 00. 

The number of noncompact components of 72 in Ml equals to the half of the 
number of its intersection points with boundary of a large rectangle {ei < / < 
e 1 _1 ,e 2 < g < e^O < £ 2 -C ei -C 1}- These intersection points correspond to the 
above asymptotic solutions and, therefore, their number equals to 4. Thus, the 
number Ni of unbounded components of 72 in is 2. □ 

2.3.4. Second application of Rolle-Khovanskii theorem. In order to estimate N 2 we 
apply Theorem l2 . 1 II again. The set T = {771 A 772 = 0} is a real algebraic curve given 
by the equation Q — 0, where 

(2.11) Qdf Adg = fg^ 3 - m A 7? 2 
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is a polynomial in (/, g). Applying Theorem 12 . 1 II again we get the following esti- 
mate. 

Proposition 2.15. 

N 2 < N 3 +N 4 , 

where (as above) N 2 is the number of points in {72 firnR^.}, -/V3 is the number of 
noncompact components ofT in M.^_ and N4 is the number of points in Tn{dQ Ary = 

Notation 2.16. For any polynomial S in two variables let NV{S) denote the 
Newton polygon of S . By ' ~i! we denote the partial order by inclusion on plane 
polygons, namely, 'A -< B' means that a polygon A lies strictly inside a polygon B. 

Lemma 2.17. The set Y has no unbounded components in R 2 . Moreover, T does 
not intersect the coordinate axes except at the origin, which is an isolated point of 

r. 

Proof. Explicit computation shows that 

(2.12) Q = 7 1 ~ 1 2 xa &fr& ~ fsQl, where Q 1 = + &' 3 ^ - 

One can easily check that = = Jjj^ = and Qi(0,0) = 0. This 

implies that the Newton polygon of the polynomial fgQi lies strictly inside the 
Newton polygon of £162 £3: 

M"P(f9Qi) = {3 < P+Q < 5, 1 < p, q < 3} < NT{Q) = {2 < p+q < 6, < p, q < 4}. 

Therefore the number of unbounded components of T in R 2 ^ coincides with the 
number of unbounded components of the zero locus of £1^2 £3 in R+, the latter 
being equal to zero. 

Another proof can be obtained by parameterizing the unbounded components of 
r near infinity and near the axes as (/ = Bt ei + g — At C2 + ...). From the shape 
of the Newton polygon of Q one can show that ei/e2 is either 0, 1 or 00. Therefore, 
B should be a root of £3, A/B should be a root of £1 or A should be a root of £2, 
respectively. Since neither of them has real roots, we conclude that T has no real 
unbounded components. 

On the coordinate axes the polynomial Q equals £i£2£3 and is therefore positive 
with the exception of the origin. The quadratic form of Q at the origin, being 
proportional to £1, is definite. Therefore the origin is an isolated zero of Q and, 
therefore, an isolated point of T. □ 

Corollary 2.18. iV 3 = 0. 

2.3.5. Estimating the number N4 of points of contact between T and r\2- These 
points are zeros of the polynomial form 

RdfAdg = fg^ 3 -dQA m . 

Thus we have to estimate the number of solutions of the system Q = R = in Mi . 

We proceed as follows. Using the Bernstein-Kushnirenko theorem we find an 
upper bound on the number of solutions of Q = R = in (C*) 2 , and then reduce 
it by the number of solutions known to be outside K 2 ^ . 

The Bernstein-Kushnirenko upper bound for the number of common zeros of 
polynomials Q and R is expressed in terms of the mixed volume of their Newton 



14 



A. GABRIELOV, D. NOVIKOV, B. SHAPIRO 



polygons. In fact, in computation of this mixed volume we replace R by its difference 
with a suitable multiple of Q: this operation does not change common zeros of Q 
and R, but significantly decreases the mixed volume of their Newton polygons. 
Simple degree count shows that the Newton polygon of R is given by 

NV{R) = {2 < p + q < 10, < p, q < 6}. 

Lemma 2.19. There exists a polynomial q = q{f,g) such that the Newton polygon 
of R = R — qQ lies strictly inside the Newton polygon of R. In other words, 

MV(R) C {3 < p + q < 9, 1 < p, q < 5}. 

Proof. Our goal is to prove that all monomials lying on the boundary of J\fV(R) 
(further called boundary monomials) are equal to monomials lying on the boundary 
of a Newton polygon of some multiple of Q. We constantly use the fact that the 
boundary monomials of a product arc equal to the boundary monomials of the 
product of boundary monomials of the factors. 

First, let us replace R by a polynomial R 2 with the same Newton polygon and the 
same monomials on its boundary, but with simpler definition. We have seen above 
^MV{SgQi)<MV{Q)=MV{ixi%i%). Denote R x df Adg = fg^ 3 -d(fgQi)A m . 
Computation of degrees shows that 

NV{Ri) C {3 < p+q < 9, 1 < p, q < 5} -< AfV(R) = {2 < p+q < 10, < p, q < 6}. 
Therefore one can disregard R\ and consider only R^ = R — Rx, where 
R 2 df Adg = M26 • d(Q - fgQ x ) A dr) 2 = const fg&Za ■ A m- 

Computing the product we see that: 



R 2 = const A 



.96 I -r^- + f&) df + /f 3 (J?- - gZ' 2 ) dg 



= const -626 { 



266 + m' 3 + 9&& - 3(1 + a)fg&&] - (1 + a) fgQx 

A simple computation using |2~Sl shows that the Newton polygon of the first product 
in the figure brackets coincides with AfV(Q) = MV(^i^2^3)- The Newton polygon of 
the second product lies strictly inside oiNV(Q), as was shown in l|2.12|l . Therefore 
it does not affect boundary monomials and can be disregarded. 

The remaining terms sum to 96661 where we denote /6C3 + 56^2 + ^66 
: 1 +a)/<76^6 by 9- Up to a non-zero constant factor, its boundary monomials are 
the same as the boundary monomials of qQ: the polynomials Q and 666 have 
proportional boundary monomials by (|2.12() . 

Using these facts we conclude that for R = R — const qQ one gets AfV(R) -< 
NT{R). ' ~ □ 

2.3.6. Bernstein- Kushnirenko theorem. Applying the well-known result of we 
know that the number of common zeros of Q and R in (C*) 2 does not exceed twice 
the mixed volume of NV{Q) and MV{R). 

Recall the definition of the mixed volume of two polygons. Let A and B be 
two planar convex polygons. It is a common knowledge that the volume of their 
Minkowsky sum XA + fiB is a homogeneous quadratic polynomial in (positive) A 
and /j,: 

Vol{\A + fiB) = Vol(A)X 2 + 2Vol(A, B)\fi + Vol{B)^ 2 . 
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By definition the mixed volume of two polygons A and B is the coefficient 
Vol(A,B). 

Setting A = fJ, = 1, one gets 

2Vol(A, B) = Vol(A + B)- Vol(A) - Vol(B). 
Lemma 2.20. There are at most 28 common zeros of Q = R = in (C*) 2 . 
Proof. Simple count gives that 2Vol(AfV(Q),AfV(R)) = 28. □ 




NV(R) 



12 3 4 

FIGURE 2. Relevant Newton polygons 



1 2 3 4 5 6 



2.3.7. Common zeros of Q and R outside R+. By Lemma l2~%l the quadratic poly- 
nomials £i, £2 and £3 have two common zeros. Denote them by (/1, <7i) and (fa, (72)- 
Evidently, they are not real (e.g. since £3 is strictly positive on R 2 ). 



Lemma 2.21. (fi,gx) and ($2,92) ore solutions of the system Q 
multiplicity at least 6 each. 



R = of 



Proof. Since (fi,gi) and (72,32) are conjugate and the system Q = R = is real, 
it suffices to consider one of these points, say (/i,.gi). 

From l|2.12|) one can immediately see that Q(fi,gi) — 0. Moreover, differentiat- 
ing Qi, one can see that = ^(fi,9i) = 0, so (A, 51) a critical point 
of Q. 

Recall that R = dQ A (gU~^ + f^) d f + ~ s€a)d»), i- e - « is thc 

product of two polynomial forms each having a simple zero at (/1 ,g\). Therefore 
this point is necessarily a critical point of R as well. 

Moreover, the 1-jet of 772 at (fi,gi) equals 



J 1? ?2 



(5 " 5i )d/ - (/ - fi)dg 



i.e. is proportional to the Euler form. Therefore, the quadratic part of R at (fi,gi) 
is proportional to the exterior product of the differential of the quadratic part of 
Q at (/i,<7i) and the Euler form, so is proportional to the quadratic form of Q 
at (/i,.9i). Thus a suitable linear combination of R and Q has both linear and 
quadratic part zero at (/i,<?i), which implies that the multiplicity of (/i,<7i) as a 
solution of the system Q = R = is at least 6. □ 
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Lemma 2.22. At least six real solutions of Q = R = lie outside K^_. 

Proof. There are exactly four points where the form 772 vanishes, exactly one in 
each real quadrant. These points are evidently solutions of the system Q = R = 0. 
Consider a connected component of the curve {Q = 0} containing such a point. It 
is a compact oval not intersecting the coordinate axes. The polynomial R vanishes 
at least once on this component, namely at this point. Therefore, R should have at 
least one another zero on this oval (counting with multiplicities), also necessarily 
lying in the same quadrant. □ 

2.3.8. Final count. The number TV of points in the intersection 71 072 PlR?_ is less or 
equal 2 + + 28— 18 = 12, where 2 is the number of unbounded components of 71 in 

; is the number of unbounded components of {Q = 0}; 28 is the Kushnirenko- 
Bernstein upper bound for the number of complex solution of R = Q = in 
(C*) 2 and 18 is the number of solutions of R = Q = outside counted with 
multiplicities. Therefore, Proposition ^ . 91 and Theorem ll.5l b are finally proved. □ 

2.3.9. Comments on Theorem \1.5\ 

1. Computer experiments indicate that, except for the two solutions of the 
system £1 = £2 = £3 = 0, all the remaining 16 solutions of the system R = Q = 
can be real. However, not all of them lie in the positive quadrant: typically the 
system R = Q = has 4 solutions in each real open quadrant. This implies 
that, provided that this statement about the root configuration could be rigorously 
proved, the best estimate obtainable by the above method would be 6, very close 
to Maxwell's conjectural bound 4. 

2. An even more important observation is that there are typically only two 
points of intersection of 72 and T — {Q = 0} lying in . In other words, the 
first application of the Rolle-Khovanskii lemma numerically seems to be exact: a 
rigorous proof that there are just two points in 71 nTHR 2 ^ would imply the original 
Maxwell conjecture. 

In fact, it is enough to prove a seemingly simpler statement that the number of 
intersections of T and 71 = — P = 0} lying m R+ 1S a t most two. This 

seems to be easier since the equation defining 71, being quadratic polynomial in g, 
can be solved explicitly. The resulting two solutions g = 51,2 (/) parameterize 71, 
and the problem reduces to the question about the number of positive zeros of a 
univariate algebraic function Q(f, <7i(/)). 

3. The fact that the polynomial R can be reduced to a smaller polynomial R 
by subtraction of a multiple of Q is a manifestation of a general yet unexplained 
phenomenon: tuples of polynomials resulting from several consecutive applications 
of the Rolle-Khovanskii theorem are very far from generic, and in every specific 
case one can usually make a reduction similar to the reduction of R to R above. 

2.4. Proof of Theorem 11.61 a. From now on we will always assume that all our 
charges are positive (the case of all negative charges follows by a global sign change) . 

2.4.1. 1-dimensional case. As a warm-up exercise we will prove Theorem II. 61 b in 
the simplest case of x-axis. 

The idea of the proof is to use the limit function 

(2.13) V 00 (x)= min ((x - Xi f + »?) = lim V-^ a {x), 
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where (xi,yi), i = I,..., I are the coordinates of the i-th charge. (We assume 
for simplicity that all yi ^ 0. The general case follows by taking the limit.) The 
function Voo(x) has at most I — 1 points of non-smoothness. Denote these points 
by 7/s. 

Lemma 2.23. Convergence V x {x) = Yim a ^ oc Va 1 ^ a {x) is valid in the C 2 -class 
on any closed interval free from jj 's. 

Proof. We assume that on such an interval p\ < (1 — rf)pi, i > 2, 77 > 0. (Here 
Pi = (x — Xi) 2 + y 2 .) Therefore, Voo{x) = p\ on this interval. The first derivative 
of Voo (x) equals 



= 2 



CiO - xi) + ^2&(Pi/pi) a 1 (x-x i ) 



i=l / 

2/a-l _ 



i=2 



(piV a (x)Y 



= 2[Ci(x - Xi) + o(l)](Ci + oil))- 2 '*- 1 - 2(x - Xl ) + o(l) = VM + o(l), 

where limQ.^00 o(l) = 0. 

Computations with the second derivative are similar, but more cumbersome. □ 

Corollary 2.24. For a sufficiently large any closed interval free from 's contains 
at most one critical point of V a (x) . 

Proof. Indeed, for any sufficiently large a the second derivative (V r Q T 1 ^ a (x))", being 
close to = 2, is positive on this interval. Therefore, V a 1 ^ a (x) is convex and can 
have at most one critical point on this interval. But the critical points of V a X ^ a (x) 
are the same as the critical points of V a (x). □ 

Lemma 2.25. For any sufficiently large a a closed interval containing some 7^ 
and free from x 4 's contains at most one critical point of V a (x) . 

Proof. Note that such an interval contains exactly one jj since 7/s are separated 
by Xi's. The required result follows from the fact that V a (x) is necessarily convex 
on any such interval. Indeed, 



(V a (x)y=a(a + l)j2<iPl a - 2 ( 
i=i ^ 



A{x- Xl ) 2 



a + 1 

Since i-ii/Oon the interval under consideration, then i s necessarily smaller 
than A{x — Xi) 2 for a large enough. Thus, (V a (x))" is positive (recall that Q > 0) 
and V a (x) itself is convex. □ 

2.4.2. Multidimensional case. We start with the discussion of the critical points 
of the limiting function Voo (x) ■ 

2.4.3. Critical points of Voo(x). The function Voo{x) is a piecewise smooth con- 
tinuous semialgebraic function. Here are the definitions of critical points of such 
functions and their Morse indices adapted to our situation. 
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Definition 2.26. A point xq is a critical point of Voc(x) if for any sufficiently 
small ball B centered at xo its subset B- = {Voo(x) < Voo(xo)} C B is either empty 
or noncontractible. 

The critical point xq is called nondegenerate if B_ is either empty or homo- 
logically equivalent to a sphere. In this case the Morse index of xq is defined as 
the dimension of this sphere plus 1. (By default, dim(0) = —1.) 

Lemma 2.27. Every effective Voronoi cell of the Voronoi diagram of a generic 
configuration of positive charges contains a unique critical point of V^x). Its 
index equals the dimension of the Voronoi cell. 

Proof. Indeed, take any effective Voronoi cell S. As above let AfS(S) denote the 
set of all nearest sites of S. By definition, S intersects the convex hull of AfS(S). 
Denote this (unique) intersection point by p(S). We claim that p(S) is the unique 
critical point of Voo(x) located on S. Indeed, the function V^x) restricted to 
J\fS (S) has a local maximum at p(S) since any sufficiently small move within AfS(S) 
brings us closer to one of the nearest sites. (Here we implicitly use the genericity 
assumptions on the configuration, i.e. that there are exactly k + 1 nearest sites for 
any Voronoi cell of codimension k and that S intersects the interior of the closure 
of J\fS(S).) On the other hand, the restriction of Voo(%) to S itself has the global 
minimum on S for similar reasons. □ 




FIGURE 3. Effective and ineffective O-dimensional Voronoi cell of Voo(x). 

Remark. Figure illustrates the above Lemma 12.271 The left picture shows the 
function Voo(x) = min(/>i, P2> Pa) where the three points are located at (1,0), 
(±-^, — |), . It is related to the left picture on Figure Q showing the corresponding 
potential V a (x) for a = 1 . In this case all the Voronoi cells of the Voronoi diagram 
are effective and one sees the local maximum inside the convex hull of these points. 
On the right picture the three points are located at (0,0), (2,0), (1, |). In this 
case the O-dimensional Voronoi cell and one of 1-dimensional Voronoi cells are 
ineffective and there is no critical point at the O-dimensional Voronoi cell. This 
picture is similarly related to the right picture on Figure ^ 

2.4.4. Proof continued. In order to settle the multidimensional case we generalize 
the previous proof using the following idea. 

Main idea for zero- dimensional Voronoi cells: Near an effective zero-dimensional 
Voronoi cell of the Voronoi diagram the union of the region where V a (x) is convex 
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(and therefore has at most one critical point) and the region where V a (x) is too 
C^-closc to V^ a (5) to have any critical points asymptotically covers a complete 
neighborhood of the Voronoi cell. More exact, we compute asymptotics of the sizes 
of the above regions, and show that the first region shrinks slower than the second 
region grows. 

The following expressions for the gradient and the Hessian form (i.e., the qua- 
dratic form defined by the matrix of the 2-nd partial derivatives) of V a (x) are crucial 
for further computations: 

/ 

WVjx) = -aY^ iP r<*-\x)V Pi {x), 
Hessy Q (S)-C = a ( a + l)^C^ 2 (S)((Vft(5),a 2 -^- T ^(5)liei| 2 ))- 

i— 1 ^ ' 

Here ' • ' denotes the evaluation of the quadratic form Hess 14 (x) at £. 

We start with the case of zero-dimensional Voronoi cells. The general case will be 
a treated as a direct product of the zero-dimensional case in the direction transversal 
to the Voronoi cell and the full-dimensional case along the Voronoi cell. 

2.4.5. Zero- dimensional Voronoi cells of Voronoi diagram. Let S be a zero-dimensional 
Voronoi cell of a Voronoi diagram of a generic configuration. We can assume that 
pi(S) = ... = p n +i(S) < Pi(S) for i > n + 1. Set 

/max i= i n+1 pi(xY 

4>{x) = log — — — - 

V imnti,. ..,„+! pi(x) 

Note that <j>(x) is everywhere positive except at the origin, and is equivalent to the 
Euclidean distance to S in a sufficiently small neighborhood of S. 

Lemma 2.28. There exists 8 > so small that in the S -neighborhood U of S the 
following conditions hold: 

(1) There exists a number e > such that for any x 6 U 

min Pi(x) > e 2e min Pi(x) > e e max pAx). 

i=n+l,n+2,...,l i=l,2,...,n+l i=l,...,n+l 

In particular, <j){x) < e inU. 

(2) The absolute value of all ratios Ck{x)/ci{x) in the unique linear dependence 
Y^i=i Ci(x)V Pi{x) = is bounded by some constant T > (by the genericity 
of configuration none of Ck{S)'s vanishes on S and, therefore, in some 
neighborhood of S as well). 

(3) If the cell S is not effective, then the closure ofU can be separated from the 
convex hull of MS{S) by a hyperplane. 

We prove that for a sufficiently large the domain U is the union of two subdo- 
mains U = U\ U U 2 such that V a (x) is convex in U\, and W a (x) ^ in U 2 - 

In what follows we denote by Ck and Kk positive constants independent of a but 
dependent on the configuration and the choice of U. 

Lemma 2.29. There exists a constant K\ independent of a such that for a suf- 
ficiently large the function V a (x) has no critical points in the domain defined by 

Ws)>^}nw. 
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Proof. First, consider the case I = n+ 1. 
The condition W a (x) = implies that 

Cj \PiJ c 3 

and, taking the logarithm of the both sides, we arrive at <p(x) < logT+log "^ c '-3( c '/ c j) _ 
Thus one can take ki — logT + logmaxj j3 '(cj/cj) in this case. 

The case I > n + 1 differs by exponentially small terms. Namely, suppose that 
Pi{x) = max pi (a;). Then 

n+l , v — a— 1 

= p? +1 W a (S) = (CiVpi - + 1> (jj Va, 

where £ = J^i=n+2 Ci(Pi/pi)~ a _1 V/9,*. One can easily see that ||£|| < Cie -6 '"" 1 " 1 ^. 
Therefore, since Vp2(x), Vp n+ i(a;) are linearly independent in U, 

KiiPi/pi)- 01 - 1 - CiCi/ci| < C 2 e- e ( Q+1 ) = o(l), 
and we get the required estimate. □ 

Lemma 2.30. There exists a constant K2 independent of a such that for all suffi- 
ciently large a the function V a (x) is convex in the domain {4>{x) < ^^rjj — — }■ 

Proof. Again, start with the case I = n + 1. The gradients V/?i,i = 1, + 1, 
span the whole M". Thus, the quadratic form J^"=i Ci(Vpi ■ £) 2 > C 3 ||£|| 2 > is 
positive definite. Therefore, one gets 

1 "+l o n+l 

-He SS y Q (x) -c = E &p7 a - 2 (vpi,o 2 — t E c^r^neii 2 > 

a(a + 1) ' a + 1 

x y i— 1 i—1 

>( max Pl )- a - 2 (YCiFtotf) - 2 (" + 1 )( mi ^i--^^)' a " lmax ^^-."+iC, ||e||2 
i=l,...,n+i V^-l / a + 1 

> (. min Pi )~ a - 2 f C 3 • (e*«)— 2 - -^U U\\ 2 - 
»=i,...,n+i y a + ly 

The last form is positive definite if 

(2.14) e -( Q+2 >^) > -Q- or, equivalent^, 0(x) < 1 °g( a + 1 )-"2 - 

a + 1 a + 2 

The case i > n + 1 differs by an exponentially small term, namely by the term 

/ 



E to 



(vp^r-^riia 2 

a + 1 



<<7 6 (. min fl i )-«- 2 e-^+ 2 )||^||- 

2=l,...,n+l 



i=n+2 

Therefore, instead of (|2.14(l we get that Hess V is positive definite provided 

e -(a+2)0(x) > + Cee -e(o ! +2) ) 

a + 1 

which gives the same estimate with a different constant. □ 

Lemma 2.31. V a {x) has at most one critical point in IA. If the cell under consid- 
eration is effective then the critical point exists and is a local minimum. If the cell 
under consideration is not effective then there is no critical point in IA. 
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Proof. Indeed, in the above notation for sufficiently large a one has 

log(a + I) — k 2 ki 
a + 2 > a + 1' 

Thus U is covered by two domains, {(f>(x) > ^j} and {<f>{x) < l£S ^^ L }- By 
Lemma \2 . 291 V a ( x ) has no critical points in the first domain. By Lemma|^|U|V' a (2;) 
is convex and has at most one critical point in the second domain. 

In the case when the considered O-dimensional Voronoi cell is effective V a (x) actu- 
ally has a local minimum located close to that Voronoi cell: the function V a (x), 
being C°-close to Voo, has a local minimum inside U. 

The last statement is a particular case of the Lemma T2.35I below. □ 

Taken together, this proves that for a sufficiently large to each effective zero- 
dimensional Voronoi cells of a generic configuration of points corresponds exactly 
one minimum of V a . 

2.4.6. Case of arbitrary codimension. Let S be any Voronoi cell of codimension 
k of the Voronoi diagram. We prove that for a generic configuration of positive 
charges and any compact K C S lying inside S there exists a sufficiently small 
neighborhood Uk independent of a containing at most one critical point of V a (x). 
Moreover, this critical point exists if and only if the cell is effective, and its Morse 
index is equal to n — k. 

Denote by L the affine subspace spanned by S. Recall that the first genericity 
assumption means that there exist exactly k+ 1 charges Ci, — >Cfe+i closest to S. 
Denote the affine subspace orthogonal to L spanned by these charges by M. 

Lemma 2.32. dimM = k. 

Proof. Indeed, a small shift of any point of S in any direction orthogonal to M 
produces a point with the same set of closest charges: distances to charges not in 
AfS(S) will still remain bigger than the distances to the charges in JVS(S), and 
the latter distances will remain equal. Therefore the shifted point still lies in S, so 
the dimension of S is at least codimM, i.e., dimM > k. The opposite inequality 
is evident since AfS(S) contains k + 1 points. □ 

If the Voronoi cell S intersects the convex hull of AfS(S) then the second gener- 
icity assumption means that any k of charges in AfS(S) do not lie on a hyperplane 
in M passing through the point L D M. 

Choosing an appropriate coordinate system we may assume that L and M in- 
tersect at the origin, i.e. are orthogonal complements of each other. Denote by xl 
and xm orthogonal projections of a vector x to linear subspaces L and M resp. i.e. 
x = xm + xl- Finally, denote the distances from x to the charges Ci, Cfc+i m 
AfS(S) by p!,...,p k +i resp. 

Let K be a compact subset of S. 

Lemma 2.33. There exists 5 > so small that in the 6 -neighborhood Uk C K™ of 
K the following conditions hold: 

(1) 3 < e <c I such that for any x G Uk one has 

min pi > e 2e min pi > e e max pi. 

i—k+2,...,l i=l,...,fe+l i=l,...,fc+l 

This is possible since K is a compact subset of an open Voronoi cell S, and 
is therefore located on some positive distance from other Voronoi cells. 
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(2) The absolute value of all ratios Ck{x)/ci(x) in the unique linear dependence 
Sj=i c k{x)^ m Pj — is bounded from above by some constant T, where 
V ' Mpj denotes the orthogonal projection of the gradient Vpj to M . Note 
that the tuple of Ck(x) is, up to proportionality, constant on S, and none of 
Ck vanishes due to the second genericity assumption. 

(3) If the cell S is not effective, then the closure ofli can be separated from the 
convex hull ofJ\fS(S) by a hyperplane. 



4> M (x) = log- 



As before, introduce the function 

' max 1=lv .. |fc+ i p 
w^t=i,...,k+i Pi 
This function is equivalent to ||xm|| near the origin: 
(2.15) C M V m (x) < \\x M \\ < CmHi(x) 

for some Cm > and all x £W. 

Lemma 2.34. a) For a certain positive K3 the function V a (x) has no critical points 
in the domain given by Uk H > K3 • e -e ( Q+1 )}. 

b) For a certain positive K4 the function V a (x) has no critical points in the 
domain given by Uk H {4>m (x) > ^j}- 

Proof, a) The idea is that outside L the gradients of pi's, i = 1, k + 1, are all 
directed away from L. The contribution of the remaining pi's, being exponentially 
small, is negligible outside an exponentially small neighborhood. 

We calculate the directional derivative of V a (x) in the direction Xl at a point 
x € Uk- 

i=l i=k+2 

Since §§^(x) = 2||5;l|| 2 for i — 1, k + 1, we conclude that the absolute value of 
the first term is at least Cg(max i= i ! ...^ + i p i )~ a ~ 1 \\xl \\ 2 - The absolute value of the 
second term is at most 

C w ( min Pi )- a - l \\x L \\ < C w ( max Pi )- a - l e-< a+ ^\\x L \\, 

i—k+2,...,l z=l,...,fc+l 

and the estimate follows. 

b) We essentially repeat the computations of Lemma \2. 291 
^ fe+i 1 

--y M v a (x) = X^r^v^Mx) + J2 Qpj^VmpAx) = 



( max 04) 

j=l,...,k+l 



i=l j=k+2 
k+1 , v -a-1 

V M p t (x) + 0(e-< a+ V) 

\maxj=i,..., k +i Pj J 



So VMV a (x) = implies, as in Lcmma r2.29l that the functions ( — 

v y 1 1 \ ^3X^=1, k+1 



Pj 



are bounded, which gives: 4>m(x) < -^fi- 



□ 



Lemma 2.35. If the cell S is not effective then for a large enough the function 
V a (x) has no critical points in Uk- 
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Proof. Let n be a direction normal to the hyperplane separating Uk from the convex 
hull of AfS(S). Then (i(Vpi,n) are all of the same sign in U (say, negative), and 
of absolute value greater than some positive constant C7. Therefore 

1 

|(W Q ,n)| = a|^CiPr _1 (V Pi ,n)|> 



> a 
a 



i=l 

C 7 (k + 1)( max p.)""- 1 - C 8 ( min p,)^ 1 

i<fe+l i>/c+l 



> 



( max PrY^ 1 ( (k + 1)C 7 - C* 8 e" e(a+1)N ) > 



,/c+l 

for a large enough. □ 

From now on we suppose that K and the convex hull of AfS(S) intersect at the 
origin, and consider the domain 

(2.16) U a = \x L < k 3 ■ e ~< a+1 \ \\x M \\ < c Uk- 

I a + 1 J 

Union of this domain and the domain described in Lemma \2 . 341 covers U by i|2.15[l . 

Moreover, <f> M {x) < in U. 

Our next goal is to study the quadratic form Hess V a {x) in the domain U a . 

Lemma 2.36. Let V^ a {x) — X)i=i CiPT"- The following holds 

a) HessV a (x) — Hess Vk, a (x) = (max i= i i ... j ( :+ i pi)~ a ~ 2 0{e~ ea ) for x G U a and 
a — > 00. 

There exist two quadratic forms A(x) and B{x) on M and L resp. such that 
Hess V k ,a(x) \-A{x)®B(x) = (max i=1 ,..., fe+1 Pj )- a - 2 0(e- ea ). 

There exist some positive constants K5, Kq such that the form A is positive definite 
and bounded from below by a k§ ■ a 2 (max.i = i^_ t k+iPi) a ' 2 , and the form B is 
negative definite and bounded from above by —kq ■ Q((niaxj = i fc+i Pi)~ a ~ 2 ■ 

Proof, a) We have to estimate from above the contribution of the distant charges. 

I RessV a (x) ■ £ - Hess V k , a (x) ■ £| < a(a + 1) £L +2 (iPi"' 2 ((V p u $ 2 - ^PiM\\ 2 ) < 
< C{ max l=1 ,..., fc+1 p % Y^e-<^ ||£|| 2 . 

b) For i = 1, k + 1 the charges Q are in M. Therefore we have 

(Vp,(x),e) 2 = (VMA(a;),6f) 2 + 2(VMP l (5),a/)(VLp l (x),a) + (VLp l (x),a) 2 = 

= {V M Pi{x),i M ) 2 + 0{e-^y url 

Here we used Vi,pj{x) = 2xl, and < K3e~ ea in W Q 
Therefore, 



fc+i „ 

a + 1 



Hess V fe , Q (x) • £ = a(a + 1) CiPp" 

i=l 

fc+1 fc+1 fc+1 
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Due to linear independence of VmPj's in U a we get, exactly as in Lemma 12.301 

i 

$>(v A /Pi,6f) 2 >C 9 ||6 



^il 2 



for some Cg > 0. Therefore, one can estimate the first term from below as 
fe+i 

a(a + l)^C^.r Q ~ 2 (VpMAf) 2 >a(a + l)C 9 (. max ft)- Q - 2 ||£A/|| 2 . 

2—1 

The second term can be estimated using 

(minCi)( max < _i_ V CpT " 1 < (max^X max p )-»-i e (°+i)*"(s). 

j'=l fe+l K + 1 J j'=l,...,fe+l 

j=l,...,fc+l 

The expression e ( Q + 1 ) < /»M(2:) j s bounded in U a by a constant e c * fK4 independent of 
a. Therefore the restriction of Hess Vk, a {x) ■ £, to M has the lower bound: 



/ L1 \/ s- a -2 ( n 2(fc + l)e c « K4 max iiW CiP l n , -ea\\ \\ e , 

a(a +1)( max p,) C* 9 — 0(e ) ) Um\ 

i=l,...,fe+l \ Ct-r l I 



2 



On the other hand, the restriction of Hess Vk, a (x) • £ to L is negative definite and 
has the upper bound: 



-2a(A;+l)minCi( max (n)-"- 1 ^]] 2 (l - 0{ e -™)) 

i— l,...,fc+l 



□ 



Corollary 2.37. For large enough a and any x € £/ Q the signature of the quadratic 
form Hess V a (x) is (k,n — k). 

Lemma 2.38. For sufficiently large a there is at most one critical point ofV a (x) 
in Uk, and its Morse index is equal to n — k. 

Proof. As was proved in Lemma 12.341 there are no critical points of V a in Uk \ U a 
for a large enough. So it is enough to prove that the mapping dV a is one-to-one in 

Take any segment / = {at — a + 1£, < t < 1} C U a . Let be the projection 
of the point dV a (a t ) on the direction £ = £m — £l (where £ = £ m + We claim 
that 7r(t) is a monotonous function of t, and, therefore, its values at ao and ai 
cannot coincide. Indeed, using Lemma \'2. 361 one can estimate 7r'(t) from below as 

n'(t) = U, ^dV a (a t ) \ = Hess V„ (f,0 > 



> ( max pi 
\=i,...,fc+i 



-a-2 



^•a 2 ||eAf|| 2 + K 6 -a||ai| 2 + 0(e- 



> 



for a sufficiently large. 

This, by convexity of U a , immediately implies the claim of the Lemma: as- 
suming that there are two critical points and joining them by a segment we get a 
contradiction. □ 

We proved in Lemma [2.351 that non-effective cells do not create critical points 
of V a for a large enough. To finish the proof of the one-to-one correspondence 
between effective cells and critical points of V a we have to show that if K contains 
a critical point of the function Voo then Uk does contain a critical point of V a ■ 
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Lemma 2.39. Assume that K contains a critical point of Voo. Then Uk contains 
a critical point of V a for a sufficiently large. 

_ j_ 

The idea of the proof is that the smooth function V a a is arbitrarily C°-close to 
Voo for large a. The topology of the sets X c = {Voo < c}(~)Uk changes as c passes the 

critical value implying the change of the topology of the sets Y c — {V a a < c} <~)Uk 

for sufficiently large a which in its turn implies presence of the critical points of 

_ j_ 

V a a , the latter coinciding with the critical points of V a . 

Let Co be the critical value of Voo at the point c G K, and let c\, ci be some 
regular values, min^ Voo < c\ < cq < C2 < maxy K Voo. We assume that c\, C2 are 
so close to Co that the interval [c± , C2] contains no critical values of Voo restricted to 
the boundary of Uk- 

_ j_ 

Let S <C (c — ci)/2. For a large enough the function V a a is at least 5/2-close 
to Voo. Thus, we have 

X c ± C Y Cl+s/2 C x ci+s c y ci+ | 5 . 

These inclusions induce homomorphisms in homology groups: 

(2.17) H4X C1 ) - H*(Y cl+5/2 ) - H4X C1+S ) - H*(Y ci + ¥ ). 

The composition of the first two homomorphisms is a homomorphism induced by 
the inclusion X Cl C X Cl+ s, which is an isomorphism. Therefore the middle homo- 
morphism in l|2.17|) is surjective. 

Similarly, the composition of the last two homomorphisms is a homomorphism 

_j_ 

induced by the inclusion Y Cl+ s/2 C Y Ci + s s , which, assuming that V a a has no 
critical points, is an isomorphism as well. Therefore the middle homomorphism in 
(I2.17|l is injective. 

Summing up, we conclude that the middle homomorphism is an isomorphism, 
and Y Cl is homologically equivalent to X Cl . 

Similarly, Y C2 is homologically equivalent to X C2 . But the sets X Cl and X C2 are 
not homologically equivalent: the first is homologically equivalent to a sphere, and 
the second is contractible. Therefore Y Cl and Y C2 are also homologically different 

implying that V a a should have a critical point in Uk- 

2.4.7. Completing the proof of Theorem, \l.fA a. Since V a has no critical points out- 
side the convex hull of the charges, it is enough to consider instead of K n an open 
ball B containing all charges. 

It is easy to see that one can cover B by neighborhoods Uk as in Lemma [2.281 
and 12.331 Indeed, start from zero-dimensional Voronoi cells Si, and choose their 
neighborhoods U according to Lemma 12.281 Then choose compacts K\ in one 
dimensional Voronoi cells in such a way that their union with these neighborhoods 
covers the intersection of the union of all one-dimensional Voronoi cells with B. 
Choose neighborhoods U K i of these compacts according to Lemma 12.331 Then 
choose compact subsets Kf of two-dimensional Voronoi cells in such a way that 
(UiCf) U (UU K i ) covers the intersection of the union of all two-dimensional Voronoi 
cells with B, and so on. At the end the union of all selected neighborhoods will 
cover B. 

Each of the neighborhoods corresponding to effective Voronoi cells will contain 
one critical point of V a , and its Morse index will be equal to the dimension of 
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the Voronoi cell. The neighborhoods of non-effective Voronoi cells will not contain 
critical points of V a . 

2.5. Proof of Theorem 11.61 b. We are looking for the critical points of the func- 
tion V a defined in a linear space N 



where 6j are now orthogonal projections on N of the positions c; of the charges Q . 
The claim is that the critical points of V a are in one-to-one correspondence with 
the effective with respect to N cells of the Voronoi diagram of {Cj}- 

We can characterize the partition of N by intersections with the cells of the 
Voronoi diagram only in terms of pi. Namely, a generalized Voronoi diagram in N 
is defined as the classical Voronoi diagram in Ml. II with pi replaced by p^. a cell 
S of a generalized Voronoi diagram is the set of all points x G N with the same 
set AfS(S) = {i|Vfc f>i(x) < pk{x)}. One can immediately see that thus defined 
generalized Voronoi diagram coincides with the intersection of the original Voronoi 
diagram with N. 

It turns out that using this notation one can get the proof of Theorem ll.GI b from 
that of Theorem 1 1.61 a by a simple replacement of M™ by N, pj by pj, the charges 
Cj by their projections on N, VV Q by VV Q , and cells of the Voronoi diagram by 
cells of the generalized Voronoi diagram. Namely, since pj is just the square of the 
distance to Cj (up to a constant), exactly the same formulae for W a and HessVo, 
hold. Since pj's are radially symmetric, the condition that N intersects the Voronoi 
diagram gencrically implies the linear independence of V/5j, which guarantees the 
second property of U and Uk , etc. 

The only difference appears for the full-dimensional strata, i.e. the case of k = 
in Lemma P2.34l and Lemma F2.36l while in Theorem ll.6l a these strata do not contain 
critical points, in Theorem l 1.61 b the full-dimensional strata corresponding to strictly 
positive pj will have a critical point. This point will be necessarily unique and is a 
local maximum by Lemma \l . 341 and Lemma P2.36I (modified as mentioned above). 

We leave it as an exercise to check that the aforementioned modifications of the 
proof of Theorem 1 1.61 a produce a correct proof of Theorem ll.Bl b. 

2.6. Proof of Theorem ll.9l Denote by a? a the number of critical points of Morse 
index j of V a (x). The standard potential V\{x) is harmonic in R 3 , and, therefore, 
has no local maxima/ minima, i.e. a\ — 0. Using the Euler characteristics one can 
easily check that a? a — Z+l = a^— for any a. Therefore, for a = 1 the total number 
of critical points a^+aj+aj equals 2a\ — l+l. By Maxwell inequalities (jl.4|l one gets 



a\ <a^ < *^,see§l. Thus, a\ + a\+a\ = 2a\-l + l <l(l-l)-l + l= (l-l) 2 , 



Remark 1. Main objects of consideration in this paper have a strong resemblance 
with the main objects in tropical algebraic geometry. Namely, the potential V a (x) 
resembles an actual algebraic hypersurface while Voo(x) resembles its tropical limit. 
Also, Voronoi diagrams are piecewise linear objects as well as tropical curves. It is 
a pure coincidence? 





exactly Maxwell's estimate. 



3. Remarks and problems 
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Remark 2. What happens in the case of charges of different signs? Note that in 
a Voronoi cell of highest dimension corresponding to a negative charge the potential 
of this charge outweighs potentials of all other charges for large a, and |V Q | _1 /" 
converges uniformly on compact subsets of this cell to V^a;). Therefore it seems 
that the function defined on the union of Voronoi cells of highest dimension as 

Voo(x) = signCi • pi(x), if pi{x) = mmpj(x) 

j 

is responsible for the critical points of V a as a — > oo. 

Remark 3. Th.eorem ll.6l is similar to the results of Varchenko and Orlik-Terao on 
the number of critical points for the product of powers of real linear forms and the 
number of open components in the complement to the corresponding arrangement 
of affine hyperplanes. Is there an appropriate result? 

Remark 4- Conjecturally the number of critical points of V a (x) is bounded 
from above by the number of effective Voronoi cells in the corresponding Voronoi 
diagram. The number of all Voronoi cells in Voronoi diagrams in R™ with I sites has 
a nice upper bound. What is the upper bound for the number of effective Voronoi 
cells? Is it the same as for all Voronoi cells? 

Remark 5. Many statements in the paper are valid if one substitutes the potential 
r~" of a unit charge located at the origin by more or less any concave function 
ip(r) of the radius in R™. To what extent the above results and conjecture can be 
generalized for ^(r)-potentials? 

Remark 6. The initial hope in settling Conjecture 11.71 was related to the fact 
that in our numerical experiments for a fixed configuration of charges the number 
of critical points of V a (x) was a nondecreasing function of a. Unfortunately this 
monotonicity turned out to be wrong in the most general formulation: the number 
of critical points of a restriction of a potential to a line is not a monotonic function 
of a. 

Example 3.1. The potential V a (x) = [(x + 30) 2 + 25] _Q + [(x + 20) 2 + 49]""+ 
[(a; + 2) 2 + 144] ~ a + [(a; - 20) 2 + 49]"" + [(a; - 30) 2 + 25]"" has three critical points 
for a — 0.1, seven critical points for a — 0.2, again three critical points for a = 0.3, 
and again seven critical points as a = 1.64, and nine critical points for a > 1.7. 

Existence of such an example for the potential itself (and not of its restriction) 
is unknown. 

4. Appendix: James C. Maxwell on points of equilibrium 

In his monumental Treatise ^U] Maxwell has foreseen the development of several 
mathematical disciplines. In the passage which we have the pleasure to present to 
the readers his arguments are that of Morse theory developed at least 50 years later. 
He uses the notions of periphractic number, or, degree of periphraxy which is the 
rank of H2 of a domain in R 3 defined as the number of interior surfaces bounding 
the domain and the notion of cyclomatic number, or, degree of cyclosis which is 
the rank Hi of a domain in R 3 defined as the number of cycles in a curve obtained 
by a homotopy retraction of the domain (none of these notions rigorously existed 
then). (For definitions of these notions see JOjj section 18.) Then he actually 
proves Theorem 11.11 of til. 1.11 usually attributed to M. Morse. Finally in Section 
[113] he makes the following claim. 
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" To determine the number of the points and lines of equilibrium, let us consider 
the surface or surfaces for which the potential is equal to C, a given quantity. Let 
us call the regions in which the potential is less than C the negative regions, and 
those in which it is greater than C the positive regions. Let Vb be the lowest and V\ 
the highest potential existing in the electric field. If we make C = Vq the negative 
region will include only one point or conductor of the lowest potential, and this is 
necessarily charged negatively. The positive region consists of the rest of the space, 
and since it surrounds the negative region it is pcriphractic. 

If we now increase the value of C, the negative region will expand, and new 
negative regions will be formed round negatively charged bodies. For every neg- 
ative region thus formed the surrounding positive region acquires one degree of 
pcriphraxy. 

As the different negative regions expand, two or more of them may meet at a 
point or a line. If n + 1 negative regions meet, the positive region loses n degrees 
of periphraxy, and the point or the line in which they meet is a point or line of 
equilibrium of the nth degree. 

When C becomes equal to V\ the positive region is reduced to the point or the 
conductor of highest potential, and has therefore lost all its periphraxy. Hence, if 
each point or line of equilibrium counts for one, two, or n, according to its degree, 
the number so made up by the points or lines now considered will be less by one 
than the number of negatively charged bodies. 

There are other points or lines of equilibrium which occur where the positive 
regions become separated from each other, and the negative region acquires pe- 
riphraxy. The number of these, reckoned according to their degrees, is less by one 
than the number of positively charged bodies. 

If we call a point or line of equilibrium positive when it is the meeting place of 
two or more positive regions, and negative when the regions which unite there are 
negative, then, if there are p bodies positively and n bodies negatively charged, the 
sum of the degrees of the positive points and lines of equilibrium will be p — 1 , and 
that of the negative ones n—1. The surface which surrounds the electrical system 
at an infinite distance from it is to be reckoned as a body whose charge is equal 
and opposite to the sum of the charges of the system. 

But, besides this definite number of points and lines of equilibrium arising from 
the junction of different regions, there may be others, of which we can only affirm 
that their number must be even. For if, as any one of the negative regions expands, 
it becomes a cyclic region, and it may acquire, by repeatedly meeting itself, any 
number of degrees of cyclosis, each of which corresponds to the point or line of 
equilibrium at which the cyclosis was established. As the negative region continues 
to expand till it fills all space, it loses every degree of cyclosis it has acquired, a 
becomes at last acyclic. Hence there is a set of points or lines of equilibrium at 
which cyclosis is lost, and these are equal in number of degrees to those at which 
it is acquired. 

If the form of the charged bodies or conductors is arbitrary, we can only assert 
that the number of these additional points or lines is even, but if they are charged 
points or spherical conductors, the number arising in this way cannot exceed (n — 
l)(n — 2) where n is the number of bodies*. 

*{I have not been able to find any place where this result is proved.}. 
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Wc finish the paper by mentioning that the last remark was added by J. J. Thom- 
son in 1891 while proofreading the third (and the last) edition of Maxwell's book. 
Adding the above numbers of obligatory and additional critical points one arrives 
at the coni ecture 1 1 . 31 which was the starting point of our paper. 
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